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The large nucleon-nucleon scattering length, and the isospin approximate symmetry, are low 
energy properties of quantum chromodynamics (QCD). These entail correlations in the binding 
energies of light nuclei, e.g., the A = 3 iso-multiplet, and Tjon’s correlation between the binding 
energy of three and four body nuclei. Using a new representation of these, we establish that they 
translate into a correlation between different short-range contributions to three body forces in chiral 
effective field theory of low-energy nuclear physics. We demonstrate that these correlations should 
be taken into account in order to avoid fine-tuning in the calibration of three body forces. We relate 
this to the role of correlations in uncertainty quantification of non-renormalizable effective field 
theories of the nuclear regime. In addition, we show that correlations can be useful in assessing the 
importance of forces induced by renormalization group (RG) transformations. We give numerical 
evidence that such RG transformations can be represented effectively by adding a constant to the 
pure three nucleon contact low energy constant ce- 
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I. INTRODUCTION 

Nuclear physics is a low-energy realization of Quan¬ 
tum Chromodynamics (QCD). However, the latter is 
non-perturbative at these low-energies, leaving this fun¬ 
damental connection hidden behind seemingly unrelated 
properties. 

The renormalization group (RG) and effective field the¬ 
ory (EFT) [1] approaches have been used for a quadran- 
scentennial to give new insights to this problem, opting 
on the separation of scales in the nuclear sector. The 
fact that the nucleon-nucleon scattering lengths are large 
compared to the nuclear force range, as indicated by the 
unnaturally small deuteron binding energy, has led to the 
development of “pionless” EFT, a renormalizable the¬ 
ory of point SU(4) nucleons interacting via contact inter¬ 
actions [2]. As only the coulomb force breaks neutron- 
proton symmetry at leading order (LO), binding energies 
of iso-multiplets are correlated, and in particular of 
and ^He [3-7]. For A = 3, pionless EFT has a limit cy¬ 
cle behavior [8], which shows the need of introducing a 
three nucleon force (3NF) already at LO. A single 3NF 
parameter naturally creates a correlation between three 
nucleon observables, such as the neutron-deutron dou¬ 
blet scattering length and the A = 3 binding energy, 
i.e., recovers the previously unexplained empirical cor¬ 
relation also known as the “Phillips line” [9-11]. Once 
such 3NF is introduced at LO, no significant cutoff vari¬ 
ation is found in the binding energy of the A = 4 system 
[10], implying that 4NF is unnecessary at this order, and 
creating a correlation between the A = 3 and A = 4 
binding energies, i.e. the “Tjon line” [10-12]. Pionless 
approach is only applicable for low momentum phenom¬ 
ena, and as a result is limited to light nuclei, where the 
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binding energy B is small enough, such that typical nu¬ 
cleon momentum Qb is much smaller than the pion mass 
Qb ~ ^ m,r {^n is the nucleon mass). 

For heavier nuclei Qb ^ and thus the theory needs 
to include pions as viable degrees of freedom. 

This is accomplished via Chiral EFT (yEFT) [1]. 
yEFT follows the notion that properties of the strong 
force are linked to the approximate chiral symmetry of 
QCD, broken at low energy to the isospin symmetry 
SU(2)i^ X SU(2)i —)• SU(2)\/. The pion is identified as 
Nambu-Goldstone boson of this symmetry breaking, and 
the pion decay constant /jr = 92.4 MeV is the order pa¬ 
rameter. This indicates that the chiral symmetry break¬ 
ing scale is A^ « 47r/^ « 1 GeV. This inherent scale 
separation, Qb ^ iri^ <C A^, invites the introduction of 
yEFT for the description of the nuclear regime. yEFT 
is constructed by integrating out chiral scale degrees of 
freedom, leaving only nucleons and pions as explicit de¬ 
grees of freedom [13]. The resulting effective Lagrangian 
retains all symmetries, particularly the approximate chi¬ 
ral symmetry, of the underlying theory. Weinberg addi¬ 
tionally showed that it can be organized in terms of a 
perturbative expansion in positive powers of QjA^. The 
coefficients of the different terms in the Lagrangian are 
called low energy constants (LECs). They capture the 
integrated out physics, and in general depend upon the 
cutoff. Currently, LECs are fixed to reproduce experi¬ 
mental data. 

Since its inception, chiral EFT, has been used to derive 
nuclear forces [14], as well as nuclear Nother currents in¬ 
duced by the global symmetries of the chiral Lagrangian 
[15]. Though the non-perturbative extension of the ex¬ 
pansion has been shown to have some weaknesses, that 
induced a debate on the correct ordering of the expansion 
[16, 17], high order yEFT at relatively high order shows 
good success in predicting nuclear properties. System¬ 
atic order-by-order improvement is needed to be demon¬ 
strated in order to assess the EFT propeties of yEFT, 
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especially for nuclei heavier than ^ = 3. In the way 
yEFT is usually practiced, i.e., Weinberg power count¬ 
ing, where the expansion is about a trivial zero momen¬ 
tum fixed point, implying naive mass dimension analysis, 
the theory is non-renormalizeable, and thus cutoff varia¬ 
tion is limited [16]. 

One of the implications is that assessing the accuracy 
of predictions, as well as quantification of uncertainties, 
which in EFTs rely sometimes on cutoff variation, are 
challenging. Moreover, in the fitting process of LECs to 
observables, EFT inherent theoretical truncation error 
is combined with experimental errors [18] . Enormous 
progress has been made in the last couple of years in 
assessing errors in predictions of EFT [19-23]. While 
applying very different approaches, all these references 
cope with the problem of how to estimate the combined 
error arising from EFT truncation and simultaneous fit 
of dozens of LECs such that the model will reproduce 
the enormous amount of data from nucleon-nucleon scat¬ 
tering. In general, order-by-order analysis is performed, 
accompanied by clever ways to propagate experimental 
error. In addition, it seems that dictating “naturalness” 
of the LECs, e.g., by using a Bayesian prior, is a very ef¬ 
ficient way to get the expected EFT order-by-order con¬ 
vergence [21, 22]. 

In the current paper we concentrate on the calibration 
process of three nucleon forces within chiral EFT. This 
poses a different challenge, due to the small number of 
3NF LECs. Up to fourth order (N^LO) there are only two 
LECs that appear exclusively in 3NF. We achieve this 
by employing known LO correlations in pionless EFT, in 
particular the Tjon’s correlation and isospin invariance 
of the nuclear forces. Not only does this allow to identify 
mistakes in the calibration of many body forces, but also 
to assess higher order contributions to observables and 
to avoid fine tuning in LECs calibrations. These make 
correlations whose origin is a different EFT to the same 
theory a powerful tool in the standing problem of quan¬ 
tifying the uncertainty in theoretical predictions of the 
EFT. Finally, by observing the evolution of such correla¬ 
tions to lower cutoffs using renormalization group trans¬ 
formations, we conjecture that the breaking of such cor¬ 
relations is a signature of neglected terms in the EFT 
potential. 


II. CALIBRATION OF xEFT FORCES 


The nuclear potential derived in yEFT is systemat¬ 
ically ordered [24-30]. For example, the leading order 
(LO) includes a pion exchange force, and singlet and 
triplet two nucleon contact terms. The next-to leading 
order (NLO) contribution includes corrections to the two 

body nuclear potential suppressed by with respect 

to the LO. suppressed terms appear at N^LO, in 


which first 3NF appear [28, 29, 31-33]. yEFT at N^LO 



FIG. 1. Contact and one-pion exchange plus contact interac¬ 
tion (a), and axial current contact (b) terms of xEFT. 


includes three different 3NF topologies: (i) N-tt-N-tt-N, 
i.e., two pion exchange; (ii) N-tt-NN, i.e., two-nucleon 
contact coupled to an additional nucleon via pion ex¬ 
change; and (iii) NNN, i.e., three nucleon contact. The 
two latter, (see Fig. 1(a)) include two new low-energy 
constants (LECs) that do not appear in the two nucleon 
or pion-nucleon sector, coined Cd and ce, respectively, 
while the long range two pion exchange term depends 
on pion-nucleon scattering parameters, that appear also 
at the two body level of N^LO xEFT, Cl, C3, C4. It is 
important to note that at N^LO no new 3NF LECs ap¬ 
pear. Up to N^LO there are about 20 LECs that appear 
in the nuclear forces. The calibration of yEFT forces 
represents a multi-variable optimization problem [19] 

In Ref. [20], the authors test a simultaneous fitting ap¬ 
proach. This is in contrast to the separate, sequential 
approach. In the latter, one notes that the LECs can 
be divided into two main sections, i.e., contact interac¬ 
tions between nucleons which represent integration out of 
physics above the chiral symmetry breaking scale, and tt- 
N LECs which are of longer wavelength. The sequential 
approach fixes these groups separately. This is due to the 
different momentum sensitivity (UV vs. lower energies), 
as well as the different physical origin, i.e., tt-N LECs 
can be calibrated using pion-nucleon scattering observ¬ 
ables only, whereas NN LECs have to include NN scat¬ 
tering data to be calibrated. This kind of separation has 
motivated bayesian approach [21, 22] to separate short- 
range and long-range contributions to the calculation er¬ 
ror. Moreover, different locality properties of the forces 
can be demonstrated by testing the convergence pattern 
when applying different regulators for these terms [34]. 

Ref. [20] has tried both. The simultaneous fixing is 
shown to have smaller order-by-order uncertainty evo¬ 
lution, and better predictivity. By carefully studying 
the optimization, the authors show major differences be¬ 
tween the approaches. In particular, they show that the 
covariance matrix of the sequential approach is nearly 
diagonal, while strong cross-correlations appear in the 
simultaneous approach. This is related to the fact that 
the experimental data of the tt-N sector is not constrain¬ 
ing well enough the theory, due to too large experimental 
errors. In addition, it is clear that one cannot separate 
the tt-N effects in the NN sector, thus the latter can be 
used also to constrain the former. However, all these 
observations can also hint to possible fine tuning. EFT 
is susceptible to fine tuning due to the fact that it is a 
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truncated expansion, and thus has an inherent error. In 
the following we focus on the role of correlations in fine 
tuning of 3NF LECs. 


Sequential fixing of 3NF LECs has been accomplished 
by the use of different observables. The triton b.e., 
together with the neutron-deuteron doublet scattering 
length [32], or with ^He b.e. [35]. Heavy nuclei properties 
were also used [36]. Lately, the fact that cd appears in 
the weak interaction operator was used to calibrate this 
3NF LEG, using the triton b.e. and decay rate as two 
constraints [37, 38]. Albeit the half-life is not as accu¬ 
rately measured as the aforementioned observables, this 
procedure has led to unprecedented accuracy of cd and 
ce fixing, a fact which we argue is related to the fact that 
the strong observables are correlated. 


The origin of the connection between the weak inter¬ 
action and 3NF is SU(2)ij x SU(2 )l symmetry, which 
is common to the chiral symmetry and the weak gaug¬ 
ing. In the strong sector, one uses the chiral symmetry 
of the yEFT Lagrangian to derive its Nbther currents 
[15]. These too are ordered according to their relative 

importance, as LO, NLO ((^a) suppressed), etc. As 
this is exactly the weak interaction gauging, these cur¬ 
rents couple to weak interaction probes. In particular, 
weak reaction rates depend, up to kinematical factors, 
upon matrix elements of these currents between the ini¬ 
tial and final nuclear states. Thus, these currents, to¬ 
gether with wave functions can be used to calculate, di¬ 
rectly from yEFT weak processes. It is instructive to 
understand the relevant structure of these currents in 
the nucleus. Up to (for currents this is N^LO), 

an axial probe interacts with a single nucleon within the 
nucleus. However, at (N^LO), two new topolo¬ 

gies appear: an interaction with a nucleon with a si¬ 
multaneous pion exchange with a different nucleon, and 
a topology in which the axial probe interacts with two 
nucleons at contact (see Fig. 1(b)). The latter is de¬ 
rived from the same term in the Lagrangian as the cd 
three nucleon force, and is thus related to cd through 
the expression [37], dn = cd + -\- 2c4) -I- 

Here, C 3 and C 4 are LECs of the dimension-two tt-N La¬ 
grangian. This relation allows the use of weak observ¬ 
ables which are sensitive to the axial current, for the cal¬ 
ibration of Cd. In particular, the decay rate of triton has 
been used to extract the axial electric J = 1 multipole 
|(^HjjE(^|j^He)| = 0.6848(11) transition matrix element 
between ^H and ^He wave functions, which are calcu¬ 
lated consistently in yEFT [37]. In the following we use 
this method, which has proven successful in many recent 
calculations (e.g., [39-41]). 


III. CORRELATIONS IN yEFT CALCULATION 
OF A = 3, 4 BOUND STATES 

The fact that yEFT and pionless EFT are EFTs of 
the same fundamental theory, means that their predic¬ 
tions should coincide at the relevant energy regimes, in 
particular Tjon’s correlation and A = 3 iso-multiplet cor¬ 
relation. The two correlations are not found in nature, 
but in computational realizations of these bound states 
using different potentials. In principle, they should be 
revealed by changing the EFT cutoff, similar to pionless 
EFT. The limited cutoff variation of yEFT potentials, 
make it impossible to check whether limit cycle behaviour 
exists in the A = 3 binding energies, and the relevance 
of the expansion about a trivial fixed point. However, 
systematically reducing the cutoff is possible using RG 
transformation [42]. Using such a method, it has been 
shown that the nuclear forces actually exhibits a Tjon 
line [35]. In particular, using two body forces only, all 
reproducing the nucleon-nucleon scattering data, leads 
to a scatter plot along a Tjon line, correlating the trinu¬ 
clei and alpha binding energies (BEs), and showing that 
neglecting three-body forces leads to at least a 20 % error 
in reproducing these nuclei BEs. In principle, a differ¬ 
ent way to generate a Tjon line is by creating different 
parameterizations of the EFT potentials at a fixed order 
and regulator, in the spirit of Ref. [19, 20]. In any case, 
yEFT leading 3NF should choose the correct BEs on this 
Tjon Line, up to higher orders. 

The Tjon line is usually plotted as a correlation line 
between ^H and "‘He BEs. In the following we show that 
for a specific choice of the two-nucleon potential, A = 3, 4 
correlations can be represented as a direct correlation 
between the three body short range LECs ce and ce- 

The BE of the '^H depends on the 3NF. We thus, for 
a specific choice of two body parameters, scan the (cd, 
Ce) plane, and search for points which reproduce the ^H 
binding energy. As seen in Fig. 2, this results in a line 
in this plane [43]. We repeat this process for ^He and 
find a different line, which follows the trend of the ^H 
line, though not overlapping. This means that while ^H 
and ^He BEs are correlated, as expected from pionless 
EFT, they cannot be described simultaneously by yEFT. 
When using EFT at a specific order v, predictions have 

inherent inaccuracy, of the order of . Thus one 

interprets the difference between triton and ^He as higher 
order effects. The difference between the lines can be 
used to asses the EFT systematic uncertainty in fixing 
these LECs. 

Indeed, this understanding is backed by studies of the 
origin of the ^H and ^He BE difference, e.g.. Ref. [3], 
where 85% of this difference is found to be a result of the 
coulomb repulsion between protons in ^He, which is de¬ 
scribed at this order of yEFT. The rest of the difference 
is found to be a result of isospin and charge symmetry 
breaking effects, specifically the proton-neutron mass dif¬ 
ference and the difference between neutron-neutron and 
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FIG. 2. Lines in {co,ce) plane which reproduce A — 3, 4 bound states, and triton strength. line is augmented to 
represent experimental triton half life uncertainty (uncertainty in BEs is smaller than the line width). The lines are plotted 
with nucleon-nucleon potentials from Ref. [29] (left pane A = 500 MeV, middle pane A = 600 MeV), and from Ref. [19] (right 
pane). 


proton-proton scattering lengths, which are of higher or¬ 
der in xEFT. 

We now repeat the same process for ^He, i.e., scan the 
plane (c_d, ce) and plot a line that reproduces ^He BE. 
In Fig. 2, we do this for three specific choices of two body 
force. For all, the difference between this line is of the 
order of the difference between the ^H-^He lines. As the 
latter is a higher order effect, we conclude that: (i) ^He 
BE is correlated to the trinuclei BE, stemming deus ex 
machina from the xEFT formalism and parameter space; 
(ii) The difference between the experimental "^He BE and 
the one predicted by a force reproducing trinuclei BE is 
a higher order effect, i.e., it is suppressed by the EFT 


IV. THE EFFECT OF RG TRANSFORMATIONS 
ON THE TJON LINE 

The use of RG transformed potentials, in which mo¬ 
menta lower than a cutoff scale A are decoupled from the 
higher ones using exact RG transformations, has become 
common in recent years [35, 42, 45]. These potentials 
do not include hard core, and thus lead to improved nu¬ 
merical convergence, and allow the use of perturbation 
theory, while exactly preserving the values of low-energy 
observables calculated with the original potential. As 
such, these potentials are ideal for many body calcula¬ 
tions. However, such a transformation induces spurious 
many body forces. In practice, one usually neglects these 
many body forces at some cluster size (usually clusters 
above A=3-4 are neglected). 

Assessing the effect of neglecting such cluster forces 
can be problematic. Here we show that prior knowledge 
regarding correlations can be used to do so. We do this by 


expansion parameter, which in this case corresponds to 
the ratio of the effective range to the scattering length. 

As an example of an observable which is not correlated, 
we plot in the (cejCe) plane the region which recovers 
the El strength of the triton beta decay. As can be 
seen in Fig.2, in the {cd, Ce) plane such a constraint is 
orthogonal to the strong BEs constraint. We augment 
the line representing the constraint by the experimental 
error. Other sources of error, such as cutoff dependence 
of the current operator [44], are not included here. The 
cutoff function and the cutoff value we use for the current 
is identical to that of the NN-tt-N 3NF diagram. 


using a viow k type of potential [42], which is constructed 
using exact RG transformation with a momentum cutoff 
A, to evolve a microscopic 2-body potential optimized to 
fit nucleon-nucleon scattering data. In this study we used 
viowk based on the 500 MeV EM potential [28]. In Fig. 3 
the correlations representation discussed in the previous 
section is given for a set of viowk evolved to differing 
momenta. One can observe that as the cutoff is reduced, 
Tjon’s correlation breaking grows. It is interesting to 
note that the correlation between the A = 3 nuclei is 
not broken, a fact that might suggest that it holds to 
higher orders. Both observations are consistent with a 
need to include three body forces induced by the RG 
transformation. 

An interesting effect we find is that Tjon’s correlation 
is broken in a way that can be fixed by adding a constant 
to Ce leg. As Ce represents the short range NNN ver¬ 
tex strength, it is reasonable that it would be affected the 
most by RG transformations, which are constructed to af- 
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FIG. 3. Lines in {cd,ce) plane which reproduce A = 3, 4 bound states, and triton strength. line is augmented to 
represent experimental triton half life uncertainty (uncertainty in BEs is smaller than the line width). Lines are plotted with 
Viowk flow of of the nucleon-nucleon potential from Ref. [28] cutoff at A = 1.6,1.8, 2.0, 2.2, 2.4, 2.6 fm“^. 


feet only short distances (similar conclusions are known 
in the NN sector, e.g., Ref. [46]). Whether a constant 
shift can be used as an effective way to avoid the techni¬ 
cal complication of including 3NF induced by RG trans¬ 
formation is a matter of future research, which should be 
accomplished by observing the behavior of different ob¬ 
servables. We note that, phenomenological evidence for 
the importance of ce for bulk nuclear properties was also 
found in Ref. [47]. 

In Fig. 3 we plot also the region in (ce, ce) plane that 


V. SUMMARY 


Pionless EFT and yEFT correspond to the same fun¬ 
damental theory, and should therefore coincide at very 
low energies. This observation allows utilizing predic¬ 
tions of one EFT as constraints to the other. In partic¬ 
ular, pionless EFT predicts at LO correlation between 
the binding energies of and ^He, and between and 
“^He. We used these correlations to give a prescription 
for assessing the quality of calibration of x^FT 3NF 
parameters {ce,Ce), for a specific choice of two-body 
nucleon-nucleon potential. Given an observable that de¬ 


recovers the strength of the triton beta decay, with 
use C 3 and C 4 from Ref. [48]. We note, that in prin¬ 
ciple one should evolve the current using the same RG 
transformation. Thus, we can give only the qualitative 
observation, that this observable is indeed useful as an 
orthogonal constraint for 3NF LEGs fixing also for viow k 
potentials. We notice that this constraint gives cd values 
stabilize for small cutoffs values. This is again consistent 
with RG philosophy of leaving long range observables in¬ 
tact. 


pends on 3NE, one plots the region in the (ce, Ce) plane 
that recovers this parameter. Using this method for the 
A = 3, 4 binding energies, one finds that they are cor¬ 
related, though the lines reproducing them are not over¬ 
lapping. Thus, they are all correlated, as expected by 
pionless EET, and the difference between them is inter¬ 
preted as higher order corrections. 

This is in fact a new representation of the correlation 
between binding energies of light nuclei, including Tjon’s 
correlation. Usually, these were demonstrated by ficti¬ 
tiously changing binding energy and showing that the 
resulting ^>"^He BEs lie on top of lines. The current pre¬ 
sentation uses physical forces, as they all reproduce real 
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world observables. In fact, had these were the only ob¬ 
servables, any point on these lines would correspond to 
a legitimate nuclear hamiltonian. As such, predictions of 
other observables using the three different lines reproduc¬ 
ing the trinuclei and alpha BEs can be used to estimate 
EFT truncation errors. The BEs are well known experi¬ 
mentally, to much higher accuracy than that of the EFT 
estimated this way. Moreover, crossings of the correla¬ 
tion lines, which represent Hamiltonians that reproduce 
accurately two of these binding energies are coincidental, 
and might lead to fine tuning. Such fine tuning can re¬ 
sult in wrong predictions and underestimation of errors. 
Contrary to the correlated BEs, we showed that the use 
of uncorrelated observable, e.g., triton /3-decay rate, is 
useful as an additional constraint. In that case we ex¬ 
pect the statistical, experimental, error rather than the 
systematical, theoretical, error to be dominant. 

We thus provide a way to disentangle systematic EFT 
errors from the experimental statistical errors in the ob¬ 
servables used to calibrate the many-body forces. More¬ 
over, in the current work we demonstrate that the naive 
combination of the these two error sources, statistical and 
systematic, can lead to fine tuning in the parameteriza¬ 
tion of the nuclear Lagrangian. Such disentanglement is 
studied in literature, e.g., using Bayesian approach, si¬ 
multaneous fitting or order-by-order study, though usu¬ 
ally in the context of the problem of calibrating many 
LECs to many data points [19-23]. Here we showed 
that in the scenario of a small number of parameters, a 
fact that usually hints to correlations, a safe calibration 
method is sequential. It will be interesting to check this 
method in similar problems, like the neutron-neutron low 


energy scattering. Ml strength, neutron-deuteron scat¬ 
tering length, or pion production [38, 49]. 

In addition, we demonstrated that breaking of these 
correlations can be used as an indication to missing terms 
in the calculation. In the case of viow k evolved potentials, 
we showed that Tjon’s correlation worsens as the cutoff 
is lowered. We thus conjectured that this difference is 
due to neglected induced many body forces in an RG 
transformation of the potential. Moreover, our results 
indicate that the majority of the effect of induced many 
body forces might be effectively represented by a constant 
(cutoff dependent) shift in the short range 3NF LEG ce- 
It also seems that there is only a moderate RG flow of 
the bare weak transition operator. In future work, where 
other observables will be studied in parallel, we intend 
to check this. A positive confirmation of these obser¬ 
vations might constitute a simple and effective solution 
to the technical problem of including currents and forces 
induced by RG transformations, as well as to better un¬ 
derstanding regularization and power counting issues. 
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